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[VII-1] 



VII. NONLINEAR OPTICS -- CLASSICAL PICTURE: 

An Extended Phenomenological model of polarization: 

As an introduction to the subject of nonlinear optical phenomena, we write, in the spirit of 
Equation [ 1-4 ], the most general form of higher order terms in the phenomenological 
electric field expansion of the polarization density (which may then be inserted in Equations 
[1-3 ])as 

P< NL, (r, = e £ J [J | \dr x dt x dr 2 dt 2 (?- r, t- t- r- r 2 t- t 2 ) E p {%, t x ) E y (? 2 , t 2 ) 

+ e o L | | | | | | d *i d h d *2 dt 2 dr 3 dt 3 Xaw ( f - f i- 1 ~ K *~ \> t- t 2 ;r- r 3 , t- 
xE p (f!,/i)E y (f 2 , t 2 ) E 8 (? 3 , t 3 ) + -. 

The wave vector and frequency dependent second and third order susceptibilities 

are then defined as 

%^ Y (kj, C0j;k 2 , u) 2 ]= \ \ \ \ d R 1 dx 1 dR 2 dx 2 exp[-z k^Rj expl+i^ xj [ VII-2a ] 

Rj R 2 T, x 2 

xexp[-i k 2 R 2 ] exp[+i© 2 x 2 ] x«p T (*i- x p ^2 * 2 ) 

and 

X^ yg [kj.cOj; k 2 , co 2 ; 63,0)3]= J J J J J J dR 1 dx l dR 2 dx 2 dR 3 dx 3 exp[-z kj ■ Rj] expf-Na^ xj 

Ri R 2 R 3 X, x 2 x 3 [ VII- 2b ] 

x exp [-z k 2 ■ R 2 ] exp[+z G) 2 x 2 ] exp [-/ k 3 ■ R 3 ] exp [+ia> 3 x 3 ] 
x%ap Y8 (Ri. x,; R 2 , x 2 ; R 3 x 3 
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Thus, we may write quite generally 



26 



P a NL '( r ' l ) = 8 <>52 Hi j^ k l d(S) l <^ k 2 d ®2 eX p[' ( k l +k 2 ) r ] eX p[-' ( C0 l +W 2)' *] 

x%a^(kp co- k 2 , co 2 ) E p (k r coJ E y (k 2 ,co 2 ) 
+e oH J J \ \ J Jc/k , t/cOj dk 2 <ico 2 dk 3 dco 3 exp[/ (k 1 + k 2 + k 3 ) fj exp[-/ (cOj +C0 2 +C0 3 ) •?] 

PY 8 \ k 2 fc, (0, co, (it 

XX aw (ki.co,;k 2 , co 2 ;k 3 , co 3 )E p (k^ co) E y (k 2 , co 2 ) E 8 (k 3 .a> 3 ) + -. 

A Simple classical model of nonlinear optical response 

A simple Lorentz-Dude model is often used in the literature as a valuable guide to the 
understanding of the frequency behavior of the nonlinear dielectric response. 27 We assume 
that the potential energy of a one-dimensional nonlinear (anharmonic) oscillator may be 
written 

V(x) =V 2 (x) + V 3 (x) + V i (x) + --- = ^Mo> 2 o x 2 + ^Ma x 3 + ^Mb / + •■• [ VII-4 ] 

(See figures on next page) 

Thus, the equation of motion of a particle moving in that potential becomes 

x+Tx+(0 2 x + ax 2 + b x 3 + ---= — E(t) [ VII- 5 ] 

M 



[ VII-3 ] 



% a p y must vanish for any material that is invariant under inversion, since both P and E are vectors, and are 

(3) 

thus odd under inversion symmetry. Note also that % a p y5 for a given material has the same transformation 

(2) (3) 

properties the elastic constants of that material. The nonzero elements of X a f,y anc ' %apy8 f° r var i° us crystal 
symmetries are compiled in Y.R. Shen, The Principles of Nonlinear Optics (Wiley, New York, 1984). 

See, for example, N. Bloembergen, Nonlinear Optics (The Advance Book Program), Addison-Wesley (1992), 
ISBN 0-201-57868-9. 
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We can analyze the response of the oscillator by expanding the displacement in powers of 
the electric field E(f) — viz. 

x[t) = x {1) { i) + x (2 \t) + x°\t) + ■■■ [ VII-6 ] 

where x ( "\t) is proportional to the nth power of the field E(f). Inserting this expression 
into Equation [ VII-5 ] and equating like powers of E(?), we obtain the following 
hierarchy of equations: 

x (i \t)+rx m (t)+(0 2 n x 0) {t) = (q/M) E(t) [VH-7a] 
x (2) {t)+rx (2 \t)+(ol x (2 \t) + a[x m {t)] 2 = [ VH-7b ] 

i (3) (f)+ri (3) (f)+(O> (3, (?) + 2ax (1) (?)jc (2) (?) + ^[x (1) (?)] 3 =0 [VH-7c] 
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In the frequency domain we see that 

-co 2 * (I) (co) - z co r;t (1) (co) +co I * (1) (co) = (q/M) E (co) [ VII-8a ] 

or * (1) (co) = (q/M) E(co) [co 2 -co 2 - z'cor] = (q/M) E(co) jF(co ,co;r) [VII-8b] 

where f(u,v\ w)= [u 2 - v 2 - iv w] . 28 Thus, knowing ^ (1) (co), we may then consider 
x m (t) a driving term in Equation [ VII-7b ] - viz. 

x {2 \t)+rx (2 \t)+(0 2 ,x (2 \t)=-a[x (l \t)] 2 [VII-9] 

Therefore, in the frequency domain 

* (2) (co) = -a ^(co ,co; I) J cko' x (l \(o) x (l \(0- (o) 1 YU ' 10a ] 

which, in view of Equation [ 1 l-8b ], becomes 
x (2) ((£>) = -a(q/M) 2 iT[co ,co;r] 



x | Jco'iT[co ,co';r] ^[co ,(co- co^r] E (co / ) E (co- co 7 ) 



.[VII-10b] 



We may treat the third order terms in a similar manner. We write Equation [ VII-7c ] as 

i (3) (?)+ri <3) (?)+co^ <3) (?)=-2a^ (1) (^ <2) (?)-^[^ (1) (?)] 3 [VII-11] 



Note that near resonance 



-si"-^"^- a i (Sw r ^ gK " m; r/2) 

where 2) is the so called complex Lorentzian. 
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In the frequency domain 

* (3) (co) = -2a ^(co ,co; T) J cko' x (1 \(o) x (2) ((0- co') . VII- 12 ] 

- b f (co ,co; r) J J d(o' d(o" x (l \(o') x (l \(o") x (1 \(0- co'- co") 

Using Equations [ VII-8b ] and [ VH-lOb ], we obtain 

x (3> ((£>) = 2a 2 (q/M) 3 ;F[co ,co;f| J J dco' dco" :T[co ,(co- co');r] ^[co ,co'; T] 

x ^fco ,co"; T] ^[cOoico'- 0)"); rl E (co- co') E(co") E(co'- co") 

ff ' [VII-13] 

-b(q/M) ^[co ,co;r] J J Jco'Jco"^[co ,co';r] J [©„,©"; I] 

x iT[co ,(co- co - 0)"); r]E(co') E(co") E (co- co - co") 

In particular, for an input 

E (t) = |Ej| cos(co 1 ?+tp : ) + |E 2 |cos(co 2 ?+(p 2 ) 

= i [Ej exp(- fcOj?) + E* exp(+z'co/)] + i [E 2 exp(-z'co 2 z) + E* exp(+zco 2 ?)] 

a. Second harmonic generation is due to the terms 

x (2) (2(D l ) =-a (q/2M) 2 ^[(H^^T] ^ 2 [(a ,(Oi,T] E^coJ [ VII-15a ] 

and 

;t (2) (2co 2 ) =-a (q/2M) 2 iF[co ,2co 2 ;f] iF 2 [co ,co 2 ;f] E 2 (co 2 ) [ VII-15b ] 

b. Sum frequency generation is due to the term 
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x (2) (a),+a) 7 ) = -a (q/2M) 2 irfco n ,((0,+(0,);r| 

vi a v» ; 27 J [ VII-16a] 

x ^[co,,,©!; T] ^[co ,co 2 ; r] E^roJ E 2 (co 2 ) 

c. Difference frequency generation is due to the term 

x (2) ((0,-a) 7 ) = -a (q/2M) 2 ^[w n ,(co,-co 7 );r] 

vi a \u i „. [ VII-16b ] 

x ^[(Oq,^; T] ^ *[co ,co 2 ; r] E^co,) E*(co 2 ) 

d. Optical rectification or dc generation is due to the terms 

x (2) {0)=-a(q/2M) 2 J [co ,0; T] |iF[co , (0{, rf ^((D^ 2 [VII-17a] 

and 

x (2) {0) =-a(q/2M) 2 ^[co ,0;f] |iT[co ,C0 2 ; if |E 2 (co 2 )f [ VII-17b ] 

e. Third harmonic generation is due to the terms 

x <?) {3(0^ = (q/M) 3 [2a 2 ^[co^G^f]- b] 



xrK,3co-r] ^K^riE^) 

and 

x (3> (3(0 2 ) = (q/M) 3 {2a 2 J[a> ,2a> 2 ;I] - b} 

x ^[a) ,3o) 2 ;r] r[(o Q ,(o 2 -,r] e^(co 2 ) 

f . Intensity dependent propagation is due to the terms 
^(cOj) = (<?/2M) 3 {2a 2 T[to ,2to{,r]-b\ 

x \7 [co , co i; if r K ,to i; T] ^(rofE^) 



[VII- 18a] 



[VII- 18b] 



[VII- 19a] 
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and 

x (3) ((d 2 ) = (q/2M) 3 {2a 2 ^[co ,2co 2 ;f] - b] 

2 [VII-19b] 
x|r[Q) ,(O 2 ;r]|Y 2 K,£0 2 ;r] |E,(© 2 )| E,(© 2 ) 

g. Raman generation (inelastic scattering) involves terms like 

x (y> (2(0 1 -0) 2 ) = 2(q/M) 3 {2a 2 ^[©^(w,-©^; r] + a 2 f[(d Q ,2(0{, Y\ - b} 

x y [©o,(2©,- co 2 );r]^[co ,co i; r] f [© ,© 2 ;T] Ef(©j) E*(© 2 ) 

Note that, according to this simple anharmonic oscillator model, Raman generation may be 
enhanced by a resonance at a frequency (©, - © 2 ) ! Also note that, for this model (see 

Equation [ VII- 13a ] above), the ratio 

Y (2) (©) 

W, \ W< \ dv \ = ~ a [ M I N V) [ VII-21 ] 

x u (g>) x (<»-<» i) x to) v ' 

is a constant independent of frequency! This observation is consistent with the famous 
empirical Miller Rule which declares that the ratio 



X^( ( °3 = tOl+t02) 

x> 3 ) x>0 x> 2 ) 



A = — J ■ — T VII-22 1 



has only a weak dispersion and is almost a constant for a wide range of materials! 

Second Harmonic Generation -- Perturbation Analysis 

We may write the nonlinear , macroscopic Maxwell equations in the form 

V 2 E(r,©) - V[ V -E(r,©)] +(© 2 /c 2 e ) I • E(r,©) = -u, © 2 P (NL) (r,©) [ VII-23 ] 



[ VII-20 ] 
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V [g ■ E(f ,©)] = -V P (NL) (f ,co) [ VII-24 ] 

Suppose that we have an input driving or pump field 

E(r,a>,)= e (1) Efa) &xp(i(k l z-(0 1 t)) + c.c. [ VII-25 ] 

Then the components of the polarization with frequency 2(0 l are given by 

P (NL) (z,co 2 = 2n) = E 2 (©J £ (2) :e (1) e (1) exp[/(2^z- 2(0$] + c.c. [ VII- 26 ] 

It is a convenience to resolve this P <NL) and the resultant second harmonic field into 
longitudinal and transverse components ~ viz. 

P (NL) (z,co 2 ) = iPj NL) (z,(D 2 ) + Pi NL) (z,(D 2 ) . [ VII-27a] 

and 

^,(d 2 )= zE„(z,G) 2 ) + E ± (z,G) 2 ) [VII-27b] 

Therefore, we may write the C0 2 or second harmonic components of Equations [ VII-23 ] 
and [ VII-24 ] 

|L E± (z,to 2 ) + ^^fel E ± (z,to 2 ) + n to 2 P ± (NL) (z,to 2 ) 
az c 8 

f ° 2 , [ VII-28a ] 

+ CO, W E||(z (02)+ ^ o W 2 p(NL) (z>(D2) z= 
I C £ J 

and — 8„(co 2 ) E ll (z,co 2 ) + P,i NL) (z,co 2 )} = [VII-28b] 

az 

To satisfy these equations two conditions must hold — viz. 
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E,(z,co 2 ) = --i-P i ; NL) (z,co 2 ) 



[ VII-29a ] 



and 



3z 2 



E ± (z,co 2 ) =-\l co 2 P ± (NL) (z,co 2 ) 



[ VII-29b ] 



where 



k 2 = — J e ±K)/ e o = 2a)iV^oe±(2cGi) 

c 



[ VII-30 ] 



We now write E ± (z,(0 2 ) = e (2) £ (z,(Q 2 ) exp ( i k 2 z) where £(z,co 2 ) is a slowly varying 
function of z — vz'z. 



A_E ± ( z ,co 2 ) = -e (2) jfc 2 £(z,co 2 )-2^ 2 A<e( z ,co 2 )-^ £(z,co 2 )|exp( ik 2 z) 
~- \k 2 2 £(z,co 2 )- 2z fc 2 ^- £(z,co 2 )|exp( i /c,z) 



[ VII-31 ] 



so that Equation [ VII-29b ] may be written 



£_£( z ,a> 2 )= — ^i o) 2 e (2) :P ± (NL) (z,a) 2 )exp(-/fc 2 z) 



3z 



2/c 



_ / M-0 <°2 Yj2/„ \ -(2) 3(2)^(1) *(1) 



[ VII-32 ] 



2L 



E 2 (co 1 ) e^:f J :e w e w exp(iMz) 



where AA= 2/q- & 2 = — i [n FH ((0,)- n SH (2a)!)] 
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If we assume that the driving field stays constant, we can directly integrate 
Equation [ VII-32 ] to obtain the exceedingly famous and important equation for 

the spatial variation of the second harmonic field — viz. 



£(2,0)2) = z 



'M-o ft>2 -(2). =(2).-(l) -(1) 
4L 



e (2) : x (2) :e [l) e w exp(i Ak z/2) ^ E 2 (cd 1 ) sinc[Mz/2] .[ VII-33 ] 



maker 29 Fringes 




-40 -20 20 

ANGULAR ROTATION IN DEGREES 



40 



LASER BEAM 



SECOND 
HARMONIC 




2=6/ COS 



Second Harmonic Generation - Coupled Wave Analysis 

When the process of second harmonic generation takes place under conditions of perfect 
phase match — i.e. k 2 = 2k x — the perturbation result breaks down if the path is 
sufficiently long. Under these circumstances the pump beam will be depleted as the second 
harmonic grows and a solution of Equation [ VII-32 ] must take into account the spatial 
variation of E 2 ((d^ . To that end we assume perfect phase matching 



P.D. Maker, R.W. Terhune, N. Nisenoff, and CM. Savage, Phys. Rev. Lett., 8, 21 (1965). 
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n (2co) 



n e (2co) 




ji/2 



and rewrite Equation [ VII-32 ] as 



|-£(z,co 2 ) = i^e (2) :I (2) :e (1) e W E^,^) 
oz /q 



[ VII-34 ] 



Of course, as the second harmonic grows Equation [ VII-3 ] tells us that a nonlinear 
polarization at the pump frequency is generated — viz. 



|<*>:e (2) e w +| (a) :e 8 e (2) 



[ VII-35 ] 



Repeating the arguments associated with Equations [ VII-28 ] through [ VII-3 1 ] we may 
write 



[ VII-36 ] 
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Combining these two equations, we find an equation governing the spatial evolution of 
E^coJ -- viz. 



eW:|«:e (2) e (1) + e (1) :| (2) :e (1) e (2) 



E(z,2(a 1 )B*(z,(ii 1 ) . [VII-37] 



Equations [ VII-34 ] and [ VII-37 ] are then the coupled differential equation which 
describe the coupling of the first and second harmonic fields. Handling all the "vectorness" 
in these two equation would obscure important issues. Thus, we consider a pair of 
somewhat less complicated equations which incorporate the essence of the problem - viz. 



3 B(l>0 j = ilb2L x w E . (l>1D j 

oz /q 



[ VII-38a ] 



|_E(z,co 1 )=%^X (2) E(z,2a) 1 )E*(z,a) 1 ) 

OZ K, 



[ VII-38b ] 



To solve these equation we first write E (z,(Oi) = E (coj fj(z) and E (z,co 2 ) = E (co,) f 2 (z) 
with the boundary conditions f^O) = 1 and f 2 (0) = 0. Thus, the coupled equations reduce 
to the dimensionless form 



[ VII-39a ] 



and 



±m=im)m 



[ VII-39b ] 



where ^= z/L e and L c '= (i cof% (2) E(co 1 )]/'^ . We next separate fj(^) and f 2 (^) into 
phase and amplitude parts as 
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U9=^(9«p[''<M9] [VII-40] 
and substitute into Equations [ VII-39a ] and [ VII-39b ] — viz. 

^u 2 {k) + iu 2 {%)^U^ = i^)^A^)-^{^]) [VII-41a] 



^ Ml (^) + / Ml (^)^ 1 (§ = / Ml (^) M2 (^)exp{-/[2^ 1 (^)-^ 2 (^)]} [VII-41b] 

Equating real and imaginary parts of these equations, we find 

^ uj& = -ufe) sin [2^(9 -<> J$] [ VII-42a ] 

11,(9 ^ 2 (9 = *?(9 cos [^i(9-*2(9] [ VH-42b ] 

j- = 49 M2 (9 sin [2^(9 — 4» 2 (9] [ VH-42C ] 

fylQ = u 2 {§ cos^)-^)] [ VII-42d ] 

Since the phase enters only in the combination \\t(Q = 2^(9 -(|) 2 (9 these four equations 
reduce to three ~ viz. 

^u 2 (§=-ufc)smy(§ [VII-43a] 

Ml (9 = 11,(9 11,(9 sin V|/(9 [ VII-43b ] 
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Combining these three equations, we obtain 



,a„«K9 fr®- 2^ A 49 --^±*M = [ VII-44a , 

which is, obviously, equivalent to 

In [ u 2 {§ ufc) cos \|/(^) ] = [ VII-44b ] 

or 

w 2 (9 wf^) cos\|/(£) = const. [ VII-44c ] 

Since w 2 (9 => as ^=> the "const' must be zero and, hence, must be 7l/2 for all 
^. Thus, the original four coupled equations now reduce to two ~ viz. 

4«2(9=«i 2 (9 [VII-45a] 

and 

^^=-^11,(9 [VII-45b] 
Combining these equations, we obtain 

4§ j- 4$ + ^(Q ^ « 2 (9 = ^ K(9 + = [ VII - 46a ] 



or 
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u?(ty + u 2 2 (ty = l [Vll-46b] 

which is an assertion of the principle of energy conservation. Taking this equation 
together with Equation [ VII-45a ] we see that 

^777 = ^ [VII-47a] 



or w 2 (^)=tanh^ [ VII-47b ] 

and d In id Q = - tanhfc) dt, => w,( Q = — ^— [ VII-47c ] 

cosh(q) 

Finally, returning to the original variables, we see that 

E(z,a> 2 )= i E(co,) tanh(z/L c ) [ VII-48a ] 
/ x E(co.) 

and E ( z ,a$ = J , T : [ VII-48b ] 

cosh(z/L c ) 
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